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HISTORY ; Detailed Current Edit wreeat 
DEC ARATION §  ; Declarative Par of odule 
DACOS = Standard Double Precision Testing DACOS 
MIHSDACOS R7 = Special DACOS routine 
MTHSDACOSD - Standard Double Precision Floating DACOSD 
MTHSDACOSD_R7 = Special DACOSD routine 
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MTHSDACOS ; Double Precision Floating Point Arc-co 16-SEP-1984 01:13:4 AX/VMS Macro V04-00 Pa 1 
1-008 ° ae ae NaF 3 UMTMRTL. SREIATHDACOS «MARS 1 ” (1) 
-TITLE MTHSDACOS ; Double Precision Floating Point Arc-cosine routine 
3; (DACOS,DACOSD) 
. IDENT /1-008/ : File: ATHDACOS.MAR Edit: JCw1008 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 

® 

¥ 

+ 

® 

THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
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: FACILITY: MATH LIBRARY 
+ 
; ABSTRACT: 


; MTHSDACOS is a function which returns the double precision floating point 
; arc-cosine in radians of its double precision floating point argument. The 
; call is standard call-by-reference. 


; MTHSDACOSD is a function which returns the double precision floating point 
3; arc-cosine in degrees of its double precision floating point argument. The 
; call is standard call-by-reference. 


. 
° 
SERRATE AAA AAAAAAAAAAAARAAAAAAAEEAEARAAAAAAAAEARERAEAAEAAAAREAEEEE 


; VERSION: 01 

; HISTORY: 

; AUTHOR: 

: Peter Yuo, 29-Jun-77: Version 01 


MODIFIED BY: 
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ion Floating Point hg 19=seF =138¢ 91: 3: $ YS Race 


-SBTTL HISTORY ; Detailed Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP=11/C ROUTINE: none 
Edit History for Version 01 of MTHSDACOS 


0-2 MTHSSERROR Res Ay | to M+ teehee 
MTH$_... chan 


0 V04- 
SRTHDACOS..MAR; 1 


Changed error Orendting "Sithenion. Put error result in RO:R1 before 
—w ange Oe eee in order to allow user modify error result. 


d MOVD #0, to CLRD. se Ls Nov-77 
ated version number and copyrt “ sftiee- 
t ange, x ABYARG to MTHSK_INVARGM setae % 
Add o PSE Cirgct tye. JBS 
Declare pan -M 


e 
Added degree entry points. RNA 22-MAR=1 
H 02-Sept-1981 


iXi>=1 


—_ 
i 
So 
o 
co 
a 


JCW 14-Jul-1982 


r+ I 1g NOV-78 


y-19 
JSB entry to THSDACOS "R7. RBG a 2t-Snt-1979 
tes Pe Py on of 1 = X*2 to avoid loss of significance for 
Corrected the floating representation of .5 from #4000 to #*x4000 
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MTHSDACOS Double Precision Floating Point Arc-co 16-SEP-1984 01:13:42 VAX/VMS M v04-00 P 
1 008 BECLARATIONS j Declarative Part of Modul moet OFF: $2 MTMRTL. SRCIMTHDACOS MAR: 1 med (3) 
78 -SBTTL DECLARATIONS 3; Declarative Part of Module 
p INCLUDE FILES: 
6 4 ; EXTERNAL SYMBOLS: 
09 6 E .DSABL GBL 
-EXTRN MTHSDSQRT_RS 
B88 8 eEXTRN DLs gd foe R7 
-EXTRN MTHSDATAND_R7 
48 90 EXTRN MTHSK_INVARGMAT 
0 91 EXTRN MTHSSSIGNAL 
$0098 
0 935 ; 
000 94 ; EQUATED SYMBOLS: 
000 95 
00004080 0000 96 $0.1.0 = *F1.0 : 1.0 
00000004 0000 97 value =4 ; value.rd.r 
0000 98 
0000 99 ; 
0000 100 ; MACROS: none 
0000 101; 
4 19 3 PSECT DECLARATIONS: 
00000000 104 -PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
4 193 ; 3; program section for math routines 
0000 107 : OWN STORAGE: none 
0000 108 ; 
000 109; 
O60 110 ; CONSTANTS: 
00 111 =; 
000 6 
4 113 D_PI_OVER_2: 
68C2 A221 OFDA 40C9 3 113 -@ORD *%*0040311, *0007732, “0121081 5 "0064302 
008 116 DPI: : 
68C2 A221 OFDA 4149 38 iN8 1 -WORD *0040511, *0007732, “0121081. *0064302 
1 119 D_90: : 
00000000 00004384 O10 131 » 18: -LONG *xX000043B84, “x0 ; 90 
00000000 00004434 0018 122 ~ LONG “%X00004434, “x0 ; 180 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: double precision floating arc-cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signals: MTHS_INVARG if !X: > 1 with reserved operand in B.-4y (copied 
to the signal “mechanism, yecter yt MCH_RO/R1 LIBS$SIG 


Associated message is: ‘‘INVAL RGUMENT™. Result is res aaa “operand -0.0 
unless a user supplied (or aa same handler changes CHFSL _MCH_RO/R1. 


" 


: 
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MTHSDACOS ; Double Precision Floating Point Arc-co 16-SEP-1984 01:13: AX/VMS Macro V04-00 P 4 Mt 
1- ATHSDACOS. = Standard Double Precision'f S-SEPcloRG $1ict ide LATUNTS eaeserMoccos.mars1 P29* (dy 1- 
~SBTTL MTHSDACOS = Standard Double Precision Floating DACOS 
3+ 
3 FUNCTIONAL DESCRIPTION: 
; DACOS = double precision floating point function 
: DACOS(X) is computed as: 
; If X = 0, then DACOS(X) = PI/2. 
8 3 If X = 1, then DACOS(X) = 0. 
; If X = <1, then DACOS(X) = PI. 
0 ; IfO<x< i, then DACOS(X) = ATAN(SQRT(1-X##2)/X), 
0 ; If 1/2.< X < 1, then DACOS(X) = ATAN(SQRT((1-X)#(14x))/X). 
0 F If -1/2 < x < 6, then DACOS(X) = ATAN(SQRT(1=X##2)/X) + PI. 
40 ; If -1 < X =< -1/2, then DACOS(X) = ATA ANCSORT CC1aK) 207140 0 + Pl. 
80 3 If 1 < iX!, error. 
00 ; CALLING SEQUENCE: 
00 3 DACOS.wd.v = MTHSDACOS(x.rd.r) 
99 : INPUT PARAMETERS: 
00000004 09 LONG = 4 3; define lLongword multiplier 
00000004 99 x = 1 * LONG S Contents of x 19 the erqusent 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


6 


3 


40FC 
6D 00000000'GF 9E 0 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


ENTRY MTHSDACOS, “M<IV, R2, R3, R4, RS, Rb, R7> 

: standard call- by-reference onry 

3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET 3 tlag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 


CONAN WWN SO ONAN EW OS OOD NAME WIN OS OOO NAME WN OOO NAUSWWANOWOONOULS aw 
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NVMINIMNNOOOOCOCOCOCOOCOCOOCOOOOOOOOSSSSSSOSOOSOOOOOOOOOOOOOOOOOOOOOOO 
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- Standard Double Precision F 
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; set handler address to jacket 


; handler 
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a 
CocecocecDCo 
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case of an error in routine 

If an atrer convert signal to user PC 
and res igna 

RO/R1 = {Xi = a@value(AP) 

call special DACOS routine 

return with result in RO/R1 


Wo 
——~ 
No 


oul 


on 
Wo 


ODooooocoo 
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ision G Las tog Poine Arc-co 1 mise} r 


= Special DAC 


JSB 
RSB 


BNEQ 
x=0 


MOVD 
RSB 


:0< ix! 


TEST_FOR.1.0: 
cw 


nn = 
owen 
= 
i 
4 


; 1 =< ix! 


GEQ_T0_1.0: 
BGTR 


BP F PAIN AIAIWIIGONINPININININININIDD 2 9 4 | 2QODTOODOOOCOCOCO OOOO OOOOOOWD@C 


WRI OOO NIA NE WN $9 OD NA UE WIN 9 OD NAME WN 0 OO NOU WAR) OO OOVNOAUE WH" OW On 
o 
= 
vw 
_—) 


routine 


MTHSDACOS_R7:: 
MTHSDACOS _R9:: 
MOvO 


3 0< ixi < 1.0 


J.13 


eee R 
MTHSDACOS_R7 


RO, R6 
TEST_FOR_1.0 


D_PI_OVER_2, RO 


#*x8000, RO 
RO, S*#SD 1.0 
GEG_T0_1.0 


.——— RO 
RO, #1, R2 
#1, RO 


R2, RO 

2$ 

RO, RO 

RO, S*#SD_1.0, RO 
MTHSDSQRT_R5 

6, RO 
MTHSDATAN_R7 
(SP)+, R6 
TEST_SIGN 


ERROR 


Pi2tiS  EMTARTL. eREIATHDACOS .MAR: 1 
~SBTTL MTHSDACOS_R7 = Special DACOS routine 

Special DACOS = used by the standard routine and direct JSB call. 

CALLING SEQUENCE: 


save anything needed in RO:R7 
OV : input in RO/R1 


; RO/R1 = PI/2 
; return PI/2 if ix! = 0 


; branch to ERROR if iXi > 1.0 


return with result in RO/R1 


special DACOS routine 
Release 1 name 

save X in R6/R7 
branch if ix! > 0 


= 
compare 


Xt 
iXi with 1.0 
branch if iX! >= 1.0 
Check for possible loss of 
arent ieamee 
R2/R5 = 1 -X 
RO/R1 = 1 + X 
RO/R1 = 1 = X*2 
Join main flow 
RO/R1 = Xe22 
seed = 1.0 = Xe22 
RO/R1 = DSQRT(1-X*2#2) 
RO/R1 = DSQRT(1-X*#2)/X 
save sign of X for sign test 
RO/R1 = DATAN(DSQRT(1-X*##2)/Xx) 
restore sign of xX 
branch to TEST_SIGN 
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MTHSDACOS ; Double Precision Floating Point Arc-co 16-SEP-1984 01:13:4 AX/VMS Macro Vv04-00 Page 7 
A bo8 ATHSDACOS_R7 = Special DACOS routine BrSep-198G 11:27:44 EMTARTE caedaribAcoS.mar:1 P22 (2) 
76 4 
re Se 
7 23 : IX! = 1.0 
a 
50 7C 7 29 CLRD RO : RO/R1 = 0 
7 
th S88 
$9 54 : Test the sign of X in order to decide if add PI to the result 
4 22 3 
7 5 TEST_SIGN: 
56 B 7 8 TSTD 6 ; test the sign of x 
06 «1 ach 5 BGEQ 10$ : branch if x > 0 
50 89 AF 60 7C¢ 60 ADDD D_FI, RO 3; add PI to RO/R1 if x < 0 
05 0080 61 10S: RSB 3: return with result in RO/R1 
SH oe 
ast 64: 1< IX!, error 
OB] 566° 
6—E DD 0081 8 ERROR: PUSHL (SP) 3; return PC from JSB routine 
7E OO°SF 9A 0083 68 MOVZBL #MTHSK_INVARGMAT, -(SP) ; condition value 
50 01 OF 79 0087 69 ASHQ #15, #T, RO 3; RO = result = reserved operand -0.0 
088 70 3 gree to signal mechanism vector 
088 71 ; (CCHFSL_MCH_RO/R1) so error handler 
088 i 3 can modify the result. 
00000000°GF 02 FB 00 7 CALLS #2, G*MTHSSSIGNAL 5 signal error and use real user's PC 
09 74 3 independent of CALL vs JS 
05 44 a RSB ; return - RO restored from CHFSL_MCH_RO/R1 
09 77 
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1 ATHSDACOSD - Standard Double Precision g- SE b= 48 4 01 43:3 ERTHRTL. SREIATHDACOS.MAR; 1 vin (3) - 


2 79 ~SBTTL MTHSDACOSD - Standard Double Precision Floating DACOSD 
38 (58) 
9 4 p++ 
4 i 3 FUNCTIONAL DESCRIPTION: 
, 5 ; DACOSD = double precision floating point function 
909 ; : DACOSD(X) is computed as: 
9 9; If X = 0, then DACOSD(X) = 90. 
9 90 ; If X = 1, then DACOSD(X) = 0. 
9 91 ; If xX = -{, then DACOSD(X) = 180. 
444 3 3 If 0 < X < 1, then DACOSD(X) = ATAN(SQRT(1=-X*##2)/Xx). 

09 93; If 1/2 < X < 1, then DACOS(X) = ATAN(SQRT((1-xX)#(14x))/X). 
009 9% ; If -1/2 < x < 6, then DACOSD(X) = ATAN(SQRT(1=X##2)/X) + 180. 

09 95; If -1 < X =< -1/2, then DACOS(X) = ATAN(SQRT((1-xX)#(14x))/X) + 180. 
009 38 : If 1< iX!, error. 

009 wy 3 

009 98 ; CALLING SEQUENCE: 

3088 99 ; 

43 s 3 DACOSD.wd.v = MTHSDACOSD(x.rd.r) 

0093 2 : INPUT PARAMETERS: 
00000004 0093 Be LONG = 4 3; define conquers multiplier 
00000004 ad b2 x = 1 * LONG ; Contents of x is the argument 

0093 09 3; IMPLICIT INPUTS: none 

009 08 ; 

it +4 3; OUTPUT PARAMETERS: 

43 \} : VALUE: double precision floating arc-cosine of the argument 

0093 13 ; IMPLICIT OUTPUTS: none 

009 15 : COMPLETION CODES: none 

009 1 ; 

44 i 3; SIDE EFFECTS: 

009 19 : Signals: MTHS_INVARG if !X! > 1 with reserved operand in RO/R1 (col80ed 

09 0 ; to the signal mechanism vector CHFSL_MCH RO/R1 y LIBSSIGNAL). 

09 1 ; Associated message is: ‘INVALID ARGUMENT”. Result is reserved operand -0.0 
43 § : unless a user supplied (or any) error handler changes CHFSL_MCH_RO/%1. 
009 4 : NOTE: This procedure disables floating point underflow, enables integer 

09 5 ; overflow. 

09 § ; 

9 Py = oe 
is 
40FC 9 0 ENTRY MTHSDACOSD, “M<IV, R2, R3, R4, RS, R6, R7> 
095 1 3 standard call-by-reference entry 
95 § 3; disable DV (and FU), enable IV 
Hy MTHSFLAG_JACKET ; flag that this is a jacket procedure in 


6D 00000000'GF 9 95 MOVAB G*MTHSSJACKET_HND, (FP) 


4 13 
Point Arc-co 1 
- Standard Double Precision 


Double Precision Floatin 


on 


74 AX/VMS Macro V04-00 P 
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3 4 put panee address to jacket 


; han 


MTHSDACOSD. R7 
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case of an error in routine 

If an a convert signal to user PC 
and resi gna 

RO/R1 = {Xi = @value(AP) 
call special DACOSD routine 
return with result in RO/R1 


MTHSDACOS ble Precision Floatin Point” A. a 16-SEP-1984 3:4 AX/VMS M v04-00 Page 1 
b08 Theda DACOSD_R RP = Special DACOSD routine g-8 SEP=1984 ate 37: 34 UMTMRIL. SR SRE RTH HDACOS.MAR; 1 : (9) 
7 tg -SBTTL MTHSDACOSD_R7 = Special DACOSD routine 
te i : Special DACOSD = used by the standard routine and direct JSB call. 
A ‘8 3; CALLING SEQUENCE: 
A 47 ; save anything needed in RO:R7 
A rk : MOV asea © 3; input in RO/Ri 
OA 49 ; JSB MTHSDACOSD_R7 
A 50 ; RSB 3 return with result in RO/R1 
OAs 380° 
pA 35 MTHSDACOSD_R7:: 3 special DACOSD routine 
56 6550 9 A 4 MOTD RO, R6 : save X in R6/R7 
06 12 OOA 55 BNEQ  D_fEST_FOR_1.0 : branch if !xX! > 0 
OOA 2§ 
OA 57: 
OA 58; x =0 
i: 
50 FF64 CF 70 OOA 61 MOVD D.90, RO : RO/R1 = 90 
05 OOAD 6¢ RSB 3: return 90 if ix! = 
QOAE 6 
OOAE 64 ; 
QOAE 65 : 0 < ix! 
A 
OAE 368 D_TEST_FOR_1.0: 
50 8000 8F AA OQOA 69 Bl cary R : RO/R1 = ix! 
08 50 71 008 70 CMPD $*#$D_1 : compare {X! with 1.0 
30 18 0086 71 BGEQ b (Gee TO_ 1 8° : branch if iX: >= 1.0 
ha 
Be 74 : 0 < iki < 1.0 
$08 a , 
50 4000 - 1 opeD AA care #04000, RO 3 Check for po es sible loss of 
52 08 50 63 O0BF 379 SUBD3 RO, #1, R2 : RORY = 1 x 
50 2 60 00C3 0 ADDD #1, RO ; RO/R1 = 1 + X 
50 64 00C6 1 MULD R2, RO : RO/R1 = 1 = xX*2 
9 11 00C9 § 8 23 : Join main flow 
30 8 64 00CB ; 1$ MULD RO, RO : RO/R1 = Xe22 
0 5 63 00C 4 SUBD RO, S*°SD_1.0, RO 3: RO/R1 = 1.0 = Xee2 
Gooog00" F ae D 5 28 JSB bie Fh ; RO/R1 = DgORT(]_xs* ) 
6 D § IVD 6, R ; RO/R1 = a re - 
6 OD 4 PUSHL 3: save sign of X for boas 
eosegeee"s 16 OB 3 JSB RINPOAT ND. R7 ; RO/R1 = DATAND cDSaRT exes »/x) 
56 3 dO E MOVL (SP)+, R6 ~ 3; restore s of X 
07 «11 OE 30 BRB D_ TEST _SIGN > branch to MTEST.SIGN 
E 92 ; 
E 3g 3 1 s¢ ix! 
2s 
E 96 D_GEQ _1T0_1.0: 
Q5 13 O0E 97 oe Rat 108 or 
FF94 31 EA 98 BRW ERROR : branch to ERROR if ‘xX! > 1.0 
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14 
gleton cet ies Petne” Arc-co 16-SEP-1984 01: 
= Special DACOSD routine 6-SEP=-1984 11: 
44 
0 
202 ¢ xt = 1.0 
208 
404 
405 10$ CLRD RO 
209 
108 
rh ; Test the sign of X in order to decide 
cite 
216 D_TEST_SIGN: 
41 TSTD 6 
414 BGEQ 10$ 
415 ADDD 0_180, RO 
218 10$: RSB 
41 
a3 
41 END 
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3; RO/R1 = 0 
if add 180 to the result 


ee the a of x 
add $80 to RO/R1 if xX <0 
return with result pi RO/R1 


MTHSDACOS 
Symbol table 


C 14 
; Double Precision Floating Point Arc-co yb ats fe 9 
66-SEP-1984 


' 


0000018 R 1 
0 10 R 1 
£Q_10_1.0 0 fe 1 
00 R 1 
OVER 2 $8 R 1 
es FOR 1.0 0 AE R 1 
TEST_SIGN 00 FR 1 
R 0000081 R 1 
GEQ_T0_1.0 be8 0074 R 1 
G = 00000004 
MTHSSJACKET_HND teeeeere xX 0] 
$$S1GNAL eereeeee = =X 8600 
MTHSDACO 00000020 RG 4 
THSDACO 0000093 RG 1 
MTHSDACOSD_R7 OOOO0AS RG 4 
MTHSDACOS_R7 00000030 RG 1 
ACOS_R9 00000030 RG 3 
MTHSDATAND_R7 eeteeeere =X 0 
ATAN_R7 teeneeee =X 00 
MTHSDSQRT_R wereeeee KX 00 
MTHSK_INVARGMAT teeeeree xX 00 
$D_1. = 00004080 
TEST_FOR_1.0 OOO0003A R 01 
EST SIGN 60000078 R 01 
= 00000004 
pew nema cme wocnane + 
; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
Os 00000000 ( 0.) 00 ¢ OO.) NOPIC USR CON’ ABS 
_MTHSCODE QOOO00F9 ¢ 249.) 01 1.) PIC USR CON” REL 
Gnccesenewesesscsesserssn$ 
; Performance indicators ; 
Page faults CPU Time Elapsed Time 
Initialization 31 0:00:00.11 8:90:00.93 
Command processing 126 0:00:00.65 0:00:04.14 
prebel sable ser 3 0: 00:01 0:00:00.47 
3 sor :00:00. :00:00. 
Pass 2 9 :00:00.97 0:00: of) 
Symbol table output 4 :00:00.03 0:00:00.11 
Psect synopsis output :00:00.01 0:00: 8°00 
Cross-reference output Bp: 8:8 9 3: 0:00.00 
Assembler run totals 35 :00:02. :00:14.14 
ye working set Limit was 900 pages. 
6252 bytes (15 pages) of virtual memory were used to buffer the intermediate code. 


There were 


pages of symbol table space allocated t 
479 source Lines were read in Pass 
1 page of virtual memory was used to define 1 macro. 


hold 26 non-local ont 


1, producing 14 object records in Pass 
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LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
LCL SHR 


EXE 


local symbols. 


RD NOWRT NOVEC LONG 


D 14 
MTHSDACOS ; Double Precision Floating Point Arc-co 16-SEP-1984 v04-00 P 1 
VAX=-11 Macro Run Statistics — , writs stay eer «any miler 9: 3}: $3 HA Leow Baran THDACOS.MAR; 1 nag 13) 


SPSS Ze sno moe ae mone nraoeneaam > 


Macro Library name Macros defined 
“S25SSDUAZB:(SYSLIBISTARLET.MLB;2 0 itt 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDACOS/GBJ=OBJ$:MTHDACOS MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) #MSRC 
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